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We discuss here a self-consistent method to calculate the properties of the cold asymmetric nu- 
clear matter. The nuclear matter is dressed with s-wave pion pairs. The nucleon-nucleon (N-N) 
f—^ I interaction is mediated by these pion pairs, uj and p mesons. The parameters of these interactions 

r— ^ , are calculated self-consistently to obtain the saturation properties like equilibrium binding energy, 

^vj ' pressure, compressibility and symmetry energy. The computed equation of state is then used in 

the Tolman- Oppenheimer-Volkoff (TOV) equation to study the mass and radius of a neutron star 
^ , containing pure neutron matter. 

^ ■ PACS numbers: 21.65.+f,21.30.Fe, 24.10.Cn,26.60.-|-c 

O' 

m : 

^ '. I. INTRODUCTION 

fj ■ The properties of nuclear matter has been an area of considerable research interest for the past few decades. Such 

^ studies are of vital importance in nuclear physics, (e.g., in the context of nucleon-nucleon (N-N) interaction, structure 

^ , and properties of finite nuclei, dynamics of heavy ion collisions), astrophysics (nucleosynthesis, structure and evolution 

'""' ■ of neutron stars [l| bordering on big-bang cosmology) and also particle physics (production and interaction among 

I [ hadrons). An obvious special case of study has been the properties of symmetric nuclear matter. However, the 

h>. ■ more general study of asymmetric nuclear matter has been receiving considerable attention of late 0, S, S| due to its 

^^ [ importance in prediction of properties of exotic nuclei, the reaction dynamics of heavy ion collisions, and properties 

t~^ . of neutron stars. 

QO One of the fundamental concerns in the study of nuclear matter is the nature of the N-N interaction. The resid- 

ual N-N interaction arising from the nucleonic substructure of quarks and gluons is basically non-perturbative in 
nature. Therefore, the general approach is to self-consistently solve the problem, albeit in different ways. The differ- 
ent approaches can be broadly classified into three general types [5|,l6|, namely, the ab initio methods, the effective 
j-IT' ■ field theory approaches and calculations based on phenomenological density functionals. The ab initio nrethods in- 
~. I elude the Brueckner-Hartree-Fock (BHF) [3, M, ISl approach, the (relativistic) Dirac-Brueckner-Hartree-Fock (DBHF) 
[13 [m, [H, [H, [ill calculations, the Green Function Monte-Carlo (GFMC) [H, [H, [l^] method using the basic N-N 
interactions given by boson exchange potentials. The other approach of this type, also known as the variational ap- 
proach, is pioneered by the Argonne Group [l^ [iSf • This method is also based on basic two-body (N-N) interactions 
C^ ■ in a non-relativistic formalism with relativistic corrections introduced at a later stage. The effective field theory (EFT) 
approaches are based on density functional theories [20, U^ like chiral perturbation theory [23, |23| ■ These calculations 
involve a few density dependent model parameters evaluated iteratively. The third type of approach, namely, the 
calculations based on phenomenological density functionals include models with effective density dep endent interac- 
tions such as Gogny or Skyrme forces [2J] and the relativistic mean field (RMF) models [2^, |2a, [27|, [281 ^m, ^^, l3l| . 
The parameters of these models are evaluated by appealing to the bulk properties of nuclear matter and properties of 
closed shell nuclei. Our work presented here belongs to this class of approaches although in the non-relativistic linrit. 
The RMF models represent the N-N interactions through the coupling of nucleons with isoscalar scalar a mesons, 
isoscalar vector u mesons, isovector vector p mesons and the photon quanta besides the self- and cross-interactions 
among these mesons [2^, [23, [Sfll- There have also been recent efforts to examine the role of isovector scalar S 
mesons 31]. Although implemented at Hartree level only, these models have been very successful in simulating the 
observed bulk properties of nuclear matter including the nuclear equation of state (EOS), mass and radii of neutron 
star as well as in explaining properties of finite nuclei [20, [2J, [23] ■ Recently, the RMF theory has been extended to 
include the quasi-particle contributions in a relativistic continuum Hartree Bogoliubov theory [32| and applied to the 
study of exotic nuclei. 
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Nuclear equations of state have also been constructed using the quark meson coupling model (QMC) ^53] where 
baryons are described as systems of non-overlapping MIT bags which interact through effective scalar and vector 
mean fields, very much in the same way as in the RMF model. The QMC model has also been applied to study the 
asymmetric nuclear matter at finite temperature [34| . 

It has been shown earlier [Sa, [3g| , that the medium and long range attraction effect provided by a mesons in RMF 
theory can be simulated by s-wave pion pairs which provide the "dressing" to the nuclear matter. Similar dressing 
of pions have also been considered to study the properties of deuteron [33] and '*He [3a|- On this basis, we start 
with a relativistic Lagrangian density with ttN interaction. The short range repulsion and the isospin asymmetry 
part of the NN interaction are parametrized by two additional terms representing the coupling of nucleons with the 
Lo and the p mesons respectively. The parameters of these interactions are then evaluated self-consistently by using 
the saturation properties like binding energy per nucleon, pressure, compressibility and the symmetry energy. The 
equation of state (EOS) of asymmetric nuclear matter is subsequently evaluated and compared with existing results 
of other independent approaches available in current literature. The EOS of pure neutron matter is then used to 
calculate the mass and radius of a neutron star. We organize the paper as follows: In Section II, we present the 
theoretical formalism of the asymmetric nuclear matter as outlined above. The results are presented and discussed in 
Section III. Finally, in the last section the concluding remarks are drawn indicating the future outlook of the model. 



II. FORMALISM 



The Lagrangian for the pion nucleon system is taken as 



£ = ■(/, {i^df, - M - Gj^Lp) ^ - 2 {dt,'Pid''ifi - in^ipupi) , 



(1) 



where ij) stands for the nucleon field with mass M, (j) = Ti4>i represents the off-mass shell isospin triplet pion field 

f ~i\ 
with mass m, Ti and 7^ being the isospin and Dirac matrices respectively, 75 = ■ p, I and G is the pion- 

nuclcon coupling constant. Repeated indices indicate summation. As shown in [35l. ISGJ] . we reduce the Lagrangian of 
equation([T]) into its non-relativistic limit and the effective Hamiltonian becomes 



H(x) = HNi^) + Wmt(x) + 7^A/(x), 

where the free nucleon part Ti.N{x.) is given by 

HAr(x)-V^t(x)e^ V^x), 

the free meson part TLm{'^) is defined as 



(2) 
(3) 



rLM{^) = \[vl + {Vvr)-{v^ 



and the ttN interaction ^35| is provided by 

7^„,t(x)=i/;t(x) 



iG 



^) + m^vl] , 



G2 



'^'^■p^"'^^ 



V^(x) 



(4) 



(5) 



In equation ([3]), ijj represents the non-relativistic two component spin-isospin quartet nucleon field and the single 
particle nucleon energy operator e^, is given by e^ — {M^ — V^)^/^. 

We expand the pion field operator ipi (x) in terms of the creation and annihilation operators of off- mass shell pions 
satisfying equal time algebra as 



'/'^(x) = 



--{ai{x.y +ai(x)), 



(pi{x) = i^l-^{ai{xy - a,;(x)) 



(6) 



with energy ujx — {rri^ ~ ^x)^ ii^ the perturbative basis. We continue to use the perturbative basis, but note that 
since we take an arbitrary number of pions in the unitary transformation \J in equation ([8]) as given later, the results 
would be nonperturbative. The two pions in eq. ([5]) provide a isoscalar scalar interaction of nucleons and thus would 



simulate the effects of cr-mesons. A pion-pair creation operator given as 

^^ = ^//(k)a.(k)ta,(-k)tdk, (7) 

is then constructed with the creation and annihilation operators in momentum space and the ansatz function /(k) 
which is to be determined later through a variational procedure. 
We then define the unitary transformation U as 

C/ = e(^'-^) (8) 

and note that U, operating on vacuum, creates an arbitrarily large number of scalar isospin singlet pairs of pions. 
The "pion dressing" of nuclear matter is then introduced through the state 

\f>^U\vac>=e^^^-^^\vac>. (9) 

Next, we define the operator U{X) with an arbitrary parameter A as U{\) — e^'^^ ^^^ and the function F{'k, A) as 
F(k, A) = t/^(A) a(k) t/(A). Differentiating i^(k. A) twice with respect to A, we have the equation 

'^^^ ^ n^FiKX). (10) 

Solving the equation (fTU|) . and identifying a(k) ~ i^(k, A = 1) we obtain 

a,(k) = [/t a,(k)(7 = (cosh /(k)) a,(k) + (sinh /(k)) a,(-k)t, (11) 

which is a Bogoliubov transformation. It can be easily checked that the operator a(k) satisfies the standard bosonic 
commutation relations: 

[a.(k), a](k')] = <5.,<5(k - k'), [a|(k), aj(k')] - [«.(k), a,(k')] = . (12) 

We then proceed to calculate the energy expectation values. We consider iV nucleons occupying a spherical volume 
of radius R such that the density p = N/{^-kR^) remains constant as (TV, R) ^ oo and we ignore the surface effects. 
We describe the system with a density operator pjv such that its matrix elements are given by |35| 

Pa/3(x,y) = Tr[pN V'/3(y)Va(x)], (13) 

and 

Tr[pNN] = f p„„(x, x)dx = N = pV. (14) 

We obtain the free nucleon energy density 

In the above equation, spin degeneracy factor 7 = 2, the index r runs over the isospin degrees of freedom n, p and 
k'f represents the Fermi momenta of the nucleons. For asymmetric nuclear matter, we define the neutron and proton 
densities p„ and pp respectively over the same spherical volume such that the nucleon density p = Pn + Pp- The Fermi 
momenta /cl and the nucleon densities are related by k'^ = (Gir'^pr/j)'^- We also define the asymmetry parameter 
y = {Pn — Pp)l P- It can be easily seen that the nucleon densities pr = |(1 ± y) for t = n,p respectively. 

Using the operator expansion of equation ([5]) , the free pion part of the Hamiltonian as given in equation ^ can be 
written as 

T-iM{^) = ai{yL)'^ LOx ai(x). (16) 

This part represents the contribution due to kinetic energy of the pions to the total energy of the system. The free 



pion kinetic energy density is given by 

hk =< /|Hm(x)|/ >^-^Jdk Luik) sinh^ /(k), 



(17) 



where a;(k) — ^/k^ + rn^. In order to calculate the interaction energy density hint in the non-relativistic limit, we 
have, from equation ([5|) using e-x — M 



Knt =< I\Tr[pN H.„t(x)]|/ >~ ^ < /I : ^,(x)^,(x) : |/ > 



Using the equations 



and pTj) . we have from equation (fT8 



i^irt.t — 



G2p 



dk /sinh2/(k) 



sinh2/(k) 



2Af V (27r)3 7 tj(k) V 2 
The pion field dependent energy density terms add up to give 

i^ra — l^k ~r flint- 

Now extremising equation (|20p with respect to /(k), we determine the ansatz function 



tanh 2/(k) = - 



G^p 



I 



2M u;2(k) 



2M 



(18) 



(19) 



(20) 



(21) 



However, we note that this ansatz function yields a divergent value for hm- This happens because we have taken the 
pions to be point like and have assumed that they can approach as near each other as they like, which is physically 
inaccurate. If we bring two pions close to each other there will be an effective force of repulsion because of their 
composite structure. We therefore replace the denominator in equation (j2ip by an additional term and rewrite the 
equation ((2T|) as 



tanh 2/(k) = - 



G^p 1 

2M ' tj^ik) + ^ + au;{k)eRl>^' 



(22) 



The introduced term in the above expression for the ansatz /(k) corresponds to a phenomenological repulsion energy 
between the pions of a "pair" given by 






3a 

(2^ 



(sinh^ /(k)) e^^-^'dk 



(23) 



where the two parameters a and Rn correspond to the strength and length scale, repectively, of the repulsion and 
are to be determined self-consistently later. Thus the pion field dependent terms of the energy density becomes 
hm — hk + hint + h5, which is then evaluated as 



3 1 



„ /G^\2 

''™" 2(2^)3V2A/J ^ 

In the above equation, the integrals Ir (t = n,p) are given by 



PnJ^n I PpJ^p 



Ir = 



UJ^ 



(^ + ae«Jfe^)i/2(^ + efiJfc^ + ^Y'^ + (u^ + e^l^^) + ^ 



(24) 



(25) 



where w = w(k). 

The short range repulsion between the nucleons is known to be mediated by the isoscalar vector u mesons. Here 
we introduce the energy of repulsion by the simple form 

K = Kp\ (26) 

where the parameter A;^ is to be fixed using the saturation properties of nuclear matter as described later. We note 



that equation (j26p can arise from a Hamiltonian density given in terms of a local potential Wfl(x) as 

Hfl(x) = ^(x)t^(x) I vni^ - ymy)H{y)dy, (27) 

where, when density is constant, we in fact have 

Ac^ == / Wfl(x)dx . 



The isospin dependent interaction is mediated by the isovector vector p mesons. We represent the contribution due 
to this interaction, in a manner similar to the w-meson energy, by the term 

hp = XpPl (28) 

where ps = (p„ — pp) and the parameter Xp is to be determined self consistently as described below. 
Thus wc finally write down the binding energy per nucleon Eb of the cold asymmetric nuclear matter: 

Eb^--M (29) 

P 

where e — {h^ -\- hf + h^ + hp) is the energy density. The expression for e contains the four model parameters a, 
i?^, Aj^ and Xp as introduced above. These parameters are then determined self-consistently through the saturation 
properties of nuclear matter. The pressure P, compressibility modulus K and the symmetry energy Esym are given 
by the standard relations: 

P - P^'-HI^ (30) 

dp 



d'{e/_p) 
dp^ 



K = V ; ; (31) 



■ld\e/p) 
The effective mass M* is given by M* = M + Vg with K, = {hint + h^)/p 



Esyrn - (2 g^, j^^^. (32) 



III. RESULTS AND DISCUSSION 

We now discuss the results obtained in our calculations and compare our results with those available in literature. 
The parameters of the model are fixed by self consistently solving eqs. (|29ll32p for the respective properties of nuclear 
matter at saturation density po — 0.15 fm^'^. The method used for solving the equations is the globally convergent 
multidimensional Secant method due to Broyden 39]. At saturation density, the values of binding energy per nucleon, 
pressure and symmetry energy are chosen to be -16 MeV, and 31 MeV respectively. The value of compressibility 
modulus K at saturation density po is chosen to be 270 MeV following a procedure described below. The pion- 
nucleon coupling strength G"^ /An, the free nucleon mass M and the pion mass m are taken to be 14.6, 940.0 MeV and 
140.0 MeV respectively. 

In order to explore the parameter space substantively, we carry out the following process: First, we manually 
tune the parameter a and using Broyden's method variationally calculate the values of Rt^ and A^^. Then, using 
these parameters in eq pip , we evaluate the compressibility K. In a similar way, next we tune Rt^ manually, while 
calculating variationally a and A^^, and subsequently the compressibility K. In both of these parameter searches, we 
find that compressibility of the nuclear matter stabilizes around the value of 270 MeV. This also shows, as expected, 
that the value of X^^ is independent of the other parameters and is around the value of 3.16 fm^. It also suggests that 
the values of a and R.,^ should be around 115 MeV and 1.06 fm respectively. 

In order to further ascertain the dependence of compressibility modulus on the parameter values, we vary the K 
value over a range 210 MeV to 280 MeV and solve variationally for a, Rt^ and A^^. It may be noted that this is 
the range of the compressibility value which is under discussion in the current literature. For K values in the range 
210 MeV to 250 MeV, the program does not converge. The solutions converge to stable parameter values only in the 
range 260 MeV to 280 MeV. The results obtained in this range for symmetric nuclear matter are presented in Table ID 
We choose K — 270 MeV and the corresponding parameter values from Table |T] for our further studies. Finally we 



generalise to asymmetric nuclear matter by fixing the symmetry energy at saturation density E^ = 31 MeV in 
eq. (|32)) and by simultaneously solving the eqs. (|29ll32p . The parameter values as shown in Table |T] remain unchanged 
in this process and we obtain Xp= 0.650 fm^. As shown later, for this set of parameter values the effective mass of 
nucleons at saturation density is found to be M* /M = 0.81. 



TABLE I: 
280 MeV. 



Parameters of the model for symmetric nuclear matter evaluated for different K values in the range 260 MeV to 
Value of saturation density pQ = 0.15 fm"''. 



K 


a 


Kit A[^ 


(MeV) (MeV) 


(fm) (fm2) 


260 


17.00 


1.424 3.098 


265 


57.28 


1.207 3.131 


270 


115.26 


1.061 3.164 


275 


205.15 


0.922 3.199 


280 


367.88 


0.756 3.237 



The binding energy per nucleon Eb as a function of the density of the system is often referred to as the nuclear 
equation of state (EOS). In the Fig. [I] we present the EOS calculated for different values of the asymmetry parameter 
y. The values y = 0.0 and 1.0 correspond to symmetric nuclear matter (SNM) and pure neutron matter (PNM) 
respectively. For SNM the binding energy, as expected, initially decreases with increase in density, reaches a minimum 
at p = po E^nd then increases. In case of PNM, the binding energy increases monotonically with increasing density. 
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FIG. 1; The binding energy per nucleon Eb as a function of nucleon density p calculated for different values of the asymmetry 
parameter y. The values y = 0.0 and 1.0 correspond to symmetric nuclear matter (SNM) and pure neutron matter (PNM) 
respectively. 



In the Figl2l we show the variation of the pressure P as a function of the energy density e. For comparison, the causal 
limit P = £ is also shown in the figure. It is evident from this curve that our EOS is consistent with the causality 
condition dP jde < 1, so that the speed of sound remains lower than the speed of light. 

Next, we present comparison of our results with the results of other groups available in the literature. In Fig.[31Ja), 
we plot Eb as a function of the nucleon density p for the symmetric nuclear matter (SNM) along with results of the 
Walecka model ^6] (long-short dashed curve) , the DBHF calculations of Li et al. with Bonn A potential (short-dashed 
curve) (data for both the models are taken from [11|) and the variational A18 -I- (5v -I- UIX* (corrected) model of 
Argonne group [l9| (long-dashed curve). While the Walecka and Bonn A models are relativistic, the variational model 
is nonrelativistic with relativistic corrections and three body correlations introduced successively. Our model produces 
an EOS softer than that of Walecka, but stiffer than the others. It is well-known that the Walecka model produces 
saturation of the nuclear matter properties correctly, though with a very high compressibility modulus oi K — 540 
MeV. Our model yields nuclear matter saturation properties correctly alongwith the compressibility oi K = 270 MeV 
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FIG. 2: The pressure P as a function of energy density e for SNM. It is evident from this curve that our EOS respects the 
causality condition dP /de < 1, so that the speed of sound remains lower than the speed of light. 
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FIG. 3: Left panehThe binding energy per nucleon Eb as a function of nucleon density p for SNM. Our results are compared 
with results of DBHF calculations with Bonn A potential [ll|], Argonne group [ly] and the Walecka model [2y]. The data for 
the Bonn A and Walecka model curves are taken from [ill ]. 
Right panel: Same as|3ja), but for PNM. 



whicli is resonably close to the empirical data. In Fig. [3I^b), we plot Eb as a function of the nucleon density p for 
PNM. Similar to the SNM case, our EOS is softer than that of Walecka model, but stiffer than those of Bonn A and 
the variational model. We use this EOS to calculate the mass and radius of a neutron star of PNM as discussed later. 
The potentials per nucleon in our model can be defined from the meson dependent energy terms of eqs. (j24p . (j26p 
and (|28p . Contribution to potential from the scalar part of the meson interaction is due to the pion condensates and 
is given by Vg — {hint + ^m)/P ^^ defined earlier. The contribution by vector mesons has two components, namely, 
due to the w and the p mesons and is given by K = Kj + V^ = (^w + hp)/p. In the Figs. IH we plot Vs and Vy as 
functions of nucleon density p calculated for SNM (Fig.[3]Ja)) and for PNM (Fig.[3]Jb)) respectively. The magnitudes 
of the potentials calculated by our model are weaker compared to those produced by DBHF calculations with Bonn 
A interaction [ll| as shown in both the panels of Fig. ID In Fig. IH^b), we show the contributions to the repulsive 
vector potential due to uj mesons (short-dashed curve), p mesons (long-dashed curve) and their combined contribution 
(long-short-dashed curve). The contribution due to p mesons rises linearly at a slow rate and has a low contribution 
at saturation density. This indicates that major contribution to the short-range repulsion part of nuclear force is from 
to meson interaction. 
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FIG. 4: Left panel: The potentials Va and VL, (as defined in the text) in SNM calculated by our model are compared with those 
of the DBHF calculations with Bonn A potential '1V\. Because of isospin symmetry, Vp (see text for definition) vanishes. Both 
the scalar (solid curve) and vector (short-dashed curve) potentials produced by our calculations are weaker in magnitude. 
Right panel: The potentials in PNM calculated by our model are compared with the Bonn A results of Li et al. 11]. The 
contributions made by the oj-meson (short-dashed curve) and p-meson (long-dashed curve) mediated interactions are distinctly 
shown for comparison. 

In Fig. [51 we present the effective mass M* /M produced by our calculations (solid curve) and those with Bonn A 
interaction (short-dashed curve) and Walecka model (long-dashed curve) (data taken from ll]). As the effective mass 
has a contribution only from the scalar potential Vg, this property is independent of the asymmetry parameter y. The 
variation of M* jM is slower in our case as compared to the Bonn A and the Walecka models. The empirical value 
for effective mass in nuclear matter derived from analysis of experimental data in the framework of non-relativistic 
shell or optical models is M* jM ~ 0.7 — 0.8. We, however, get a shghtly higher effective mass of M* jM = 0.81 
at saturation density. This is due to lower contribution (in absolute terms) of the scalar interaction in the medium as 
shown in Fig. [5l 
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FIG. 5: The effective mass of nucleons (as defined in text) calculated by our model (solid curve) is compared with those of 
Bonn A (short-dashed) and Walecka model (long-dashed) curves (data taken from [llj]). At saturation density the M* jM value 
produced by our model is 0.81. A weaker potential Va produces a higher M* jM value in our calculations. 



Knowledge of density dependence of symmetry energy is expected to play a key role in understanding the structure 
and properties of neutron-rich nuclei and neutron stars at densities above and below the saturation density. Therefore 
this problem has been receiving considerable attention of late. Several theoretical and experimental investigations 
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FIG. 6: Symmetry energy Esym calculated from the EOS(as in Eq. I32p (solid blue line) is plotted along with results of other 
groups. The data for experimental points and the results of the antisymmetrized molecular dynamics (AMD) simulations with 
Gogny-AS and Gogny interactions are taken from Shetty et al ^l|], DBHF (Bonn A) results are taken from [231, RMF (NLp) 
data are frornJsB, the variational (A18+(?v+ UIX*) model of Akmal et al. il!l| results are from '^, DBHF {aujp) model of 
Huber et al. [42. k3| data are from [23], the Skyrme amd NL3 results are from _5!|. Our result shows consistency with those of 
other groups and corroborates the "stiff" dependence of Esym as advocated by Shetty et al. 1411] . 



addressing this problem have been reported (0, [^, |4l| and references therein). While the results of independent 
studies show reasonable consistency at sub-saturation densities p < po, they are at wide variance with each other at 
supra-saturation densities p > pq. This wide variation has given rise to the so-called classification of "soft" and "stiff" 
dependence of symmetry energy on density [40|, |4l| . 

Fig. [6] shows a representation of the spectrum of such results alongwith the results of the present work (solid blue 
curve) . While the Gogny and Skyrme forces (dark rib-dotted and dotted curves respectively with data taken from 
[SjEil) produce "soft" dependence on one end, the NL3 force (dot-dashed curve with data taken from Q) produces a 
very "stiff" dependence on the other end. The analysis of experimental and simulation studies of intermediate energy 
heavy-ion reactions as reported by Shetty et al. [4l| (red tri ang les and long-short-dashed red curve repectively) , results 
of DBHF calculations of Li et al. and Huber et al. [ill, [2^, E, H| (rib-dashed and magenta ribbed curve), variational 
model [J, [l0| (short-dashed curve) , RMF calculations with nonlinear Walecka model including p mesons by Liu et 
aL|3]| (long-dashed green curve) as shown in Fig. [6] suggest "stiff" dependence with various degrees of stiffness. The 
experimental results (represented by the red triangles with data taken from Shetty et al. [4l|) are derived from the 
isoscaling parameter a which, in turn, is obtained from relative isotopic yields due to multifragmentation of excited 
nuclei produced by bombarding beams of ^®Fe and ^^Ni on ^*Fe and ^*Ni targets. Shetty et al. have shown that the 
results of multifragmentation simulation studies carried out with Antisymmetrized Molecular Dynamics (AMD) model 
using Gogny-AS interaction and Statistical Multifragmentation Model (SMM) are consistent with the above-mentioned 
experimental results and suggest (as shown by the red long-short-dashed curve) a moderately stiff dependence of the 
symmetry energy on density. Our results (represented by the solid blue curve) calculated using eqn. (|32p are consistent 
with these results at subsaturation densities but are stiffer at supra-saturation densities. In Fig[6l the curve due to 
Huber et al. [43,|43] (with data taken from [2^) correspond to their DBHF 'HD' model calculations which involves 
only the a, lu and p mesons. Similarly the long-dashed green curve due to Liu et al. |31j is from the basic non-linear 
Walecka model with cr, to and p mesons. Our formalism is the closest to these two models with the exception that in 
our model the effect of a mesons is simulated by the tt meson condensates. It is also noteworthy that our results are 
consistent with these results for densities upto 2po- 

The wide variation of density dependence of symmetry energy at supra-saturation densities has given rise to the 
need of constraining it. As discussed by Shetty et al [4l| . a general functional form Esym = ^symip/Po)'^ h^-s 
emerged. Studies by various groups have produced the fits with -E° „ ~ 31 — 33 MeV and 7 ~ 0.55 — 1.05. A similar 
parametrization of the Egym produced by our EOS with E'^ — 31 MeV yields the exponent parameter 7 = 0.85. 
In Fig. [3 the symmetry energy calculated directly from the EOS by eqn. (|5^ (solid curve) and by the fit (dashed 
curve) are presented. 

We next use the equation of state for PNM derived by our model in the Tolman-Oppenheimer-Volkoff (TOV) 
equation to calculate the mass and radius of a PNM neutron star. In Fig. [5] we show the neutron star mass as a 
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FIG. 7: Density dependence of the symmetry energy Esym 
the blue sohd hne of Fig. |6]) and from parametric fit Eaym 
plotted. 



as calculated from the EOS using eqn. 



(solid curve, same as 



E%^{plp^r with E'l 



31 MeV and 7 = 0.85 (dashed) are 



function of its radius as obtained from the TOV equation. The mass and radius of the star are found to be 2.25 Mq 
and 11.7 km respectively. In the Fig. [5] we show the neutron star mass as a function of its central energy density £c- 
We observe that as Sc increases the mass of the star increases first but eventually reaches the Chandrasekher limit. 
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FIG. 8: Mass of a PNM neutron star (in units of solar mass Mq) produced by our EOS as a function of radius. Maximum 
mass of the star turns out to be 2.25 Mq with a radius of 11.7 km. 



IV. CONCLUSION 



In this work we have presented a quantum mechanical nonperturbative formalism to study cold asymmetric nuclear 
matter using a variational method. The system is assumed to be a collection of nucleons interacting via exchange of tt 
pairs, a; and p mesons. The equation of state (EOS) for different values of asymmetry parameter is derived from the 
dynamics of the interacting system in a self-consistent manner. This formalism yields results similar to those of the 
ah initio DBHF models, variational models and the RMF models without invoking the u mesons. The compressibility 
modulus and effective mass are found to be X = 270 MeV and M* jM = 0.81 respectively. The symmetry energy 
calculated from the EOS suggests a "stiff" dependence at supra-saturation densities and corroborates the recent 
arguments of Shetty et al. [41[. A parametrization of the density dependence of symmetry energy of the form 
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FIG. 9: Mass of a PNM neutron star as a function of the central energy density £c is plotted. The mass of the star converges 
to Chandrasekhar limit. 
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^'iymip/PoV with the symmetry energy E'^^^ at saturation density being 31 MeV produces 7 — 0.85. The 



EOS of pure neutron matter (PNM) derived by the formalism yields the mass and radius of a PNM neutron star to 
be 2.25 Mq and 11.7 km respectively. 

Besides the aesthetic appeal, we note that the present formalism may also have experimental consequences. For 
example, the off shell 7r"'"7r~-pair in the pion dressing may annihilate to hard photons with a probability in excess of 



what one may expect otherwise, 
is currently under progress. 



Extension of the formalism to study asymmetric nuclear matter at finite temperature 
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